Unconventional superconducting gap via spin fluctuations in iron-vacancy ordered 
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Based on an effective f2-orbital tight-binding model, we examine the superconducting states 
induced by the antiferromagnetic fluctuations for iron- vacancy-ordered Aj;Fe2-a;Se2. It is shown 
that due to the broken reflection symmetry induced by the iron vacancies, new superconducting 
states with Cah symmetry emerge. In particular, we show that in the C4,h symmetry, symmetric 
axes of the pairing momenta do not need to coincide with axes of the unit cell. As a result , in addition 
to the magnitude of the pairing gap, the relative orientation of the pairing wave function to the 
lattice forms another degree of freedom for characterizing the superconducting gap and can further 
help in gaining the condensation energy. Nonetheless, similar to other iron-based superconductors, 
the singlet ground state is still dominated by s -wave or d-wave, which are nearly degenerate 
with anisotropic gaps. Furthermore, s-wave and d-wave superconducting states are separated by a 
quantum critical point controlled by the Hund's rule coupling .Jh- 

PACS numbers: 74.70.Xa, 74.20.Mn, 74.20.Rp 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity in cuprates, it has been a key interest to find the 
mechanism that causes the high Tc and its unconven- 
tionahty. After more than 20 years of intensive investi- 
gations, however, the origin of high- Tc superconductiv- 
ity in cuprates is still unsolved. The discovery of rel- 
atively high Tc in iron-based superconductors^'^ opens 
a different route for high-temperature superconductivity 
and thus provides a unique opportunity to re-clarify the 
physics of the unconventional superconductivity. One of 
the features for the unconventional superconductivity is 
the appearance of the sign-switched gap function. The 
presence of disconnected Fermi surface (FS) sheets in 
iron-based superconductors is an ingenious condition for 
the existence of the sign-switched order parameters such 
as s±-wave and d-wave"^. It is widely believed that the 
superconducting (SC) pairing is driven by the inter-FS- 
sheet scattering or equivalently by the antiferromagnetic 
(AFM) fluctuations^"— as the inter-FS-sheet scattering 
also leads to a strong AFM correlation. 

The recent discovery of alkaline-intercalated iron- 
selenide superconductor AyFe2-a;Se2 with above 
SOKZt initiates an interesting branch for exploring su- 
perconductivity in iron-based superconductors. The 
chemically-stable structure A0.8Fe1.6Se2 exhibits a 
•\/5 X \/b iron-vacancy order and block-checkerboard 
antiferromagnetism^"— . By varying contents of iron or 
alkaline, the AFM state becomes unstable and the SC 
state appearsiiii^. More recently, it is indicated that 
the antiferromagnetic phase and the superconducting 
phase may be separated in different space regions^"— . 
In particular, a scanning tunneling microscopy (STM) 
measuremenlr^ shows that the sample would be spon- 
taneously separated into the superconducting KFe2Se2 
and the insulating Ka;Fei.6Se2. Nevertheless, it still 



needs further clarification to see if the insulating re- 
gion can be altered into a metal by chemical doping or 
by applying pressure. First-principles calculations in- 
dicate that by applying pressure to samples with the 
%/5 X -v/S iron- vacancy order, the system goes through 
two successive magnetic transitions, from the semicon- 
ducting block-checkerboard AFM phase to the metal- 
lic stripe (collinear) AFM phase and then to a metallic 
non-magnetic phase^S,. These results agree with trans- 
port measurements and indicate that a semiconductor-to- 
metal transition can be induced by the pressure^i. The 
issue of whether the iron-vacancy-ordered AyFe2-a;Se2 
supports superconductivity, however, remains unsettled 
as it is shown that as the system became more metallic 
under pressure, Tc gets smaller, and whether it suggests 
the necessity of strong correlation for the emergence of 
superconductivity is an open question. 

On the other hand, nuclear magnetic resonance (NMR) 
experiments find singlet superconductivity and no pro- 
nounced spin fluctuations near T^S"^. The spin-lattice 
relaxation rate 1/ Ti T shows the absence of the Hebel- 
Slichter coherence peak and a power-law behavior below 
Tc, indicating an unconventional SC gap which is very 
likely gapless. Although salient spin fluctuations are not 
present in NMR, the possibility of spin-fluctuation-driven 
superconductivity should not be excluded since the infor- 
mation of spin fluctuations would be hidden in signals^S. 
A recent Raman measurement on the two-magnon scat- 
tering shows that the scattering rate grows as the temper- 
ature approaches Tc and undergoes a sudden drop when 
the system enters the SC phase^^. It thus supports that 
superconductivity arises from magnetic fluctuations. 

Theoretically, previous works^^"— on superconductiv- 
ity in Aj^Fe2-a;Se2 were based on the band structure 
of KFe2Se2r^i^. However, the iron- vacancy order has 
been proved to induce a big change in the FS shap o^°'^^ , 
though it is not confirmed so far by the angle-resolved 



photoemission spectroscopy (ARPES)^"— . Therefore, if 
superconductivity could happen in the iron-vacancy or- 
dering phase, it is crucial to examine the pairing symme- 
try and pairing mechanism. 

In this work, we investigate superconducting insta- 
bility in iron- vacancy-ordered AyFe2-xSe2. Previously, 
based on the 12-orbital tight-binding model that fits the 
band structure of Kj,Fe2-xSe2, we have studied the mag- 
netic phase in the generalized Hubbard model and suc- 
ceeded in explaining the block-checkerboard AFM insta- 
bility from the Stoner's theory^. Here based on the same 
tight-binding model and starting from a non-magnetic 
metallic phase, we examine the superconductivity by the 
fluctuation-exchange (FLEX) approacb^2ii2, . For sin- 
glet superconductivity, we examine the pairing coupling 
constants for s-wave and d-wave states from the effec- 
tive scattering matrix. We find that similar to other 
iron-pnictide superconductors^, s-wave and rf-wave are 
closely degenerate. Furthermore, the gap functions in 
both symmetries are highly anisotropic and have nodes. 
In addition, since the iron- vacancy order lowers the sym- 
metry of AyFe2-xSe2 to the group C4h, the gap func- 
tions do not need to be refiection symmetric and this 
implies that there are degrees of freedom in defining 
origins of angles for the pairing momenta. These ex- 
tra degrees of freedom for the gap functions can further 
increase the gain of the condensation energy. In addi- 
tion, we find that the pairing symmetry is controllable 
by tuning the Hund's rule coupling Jh and the critical 
value at about Jh — 0-2U. The Hund's rule coupling 
manipulates the phase transition or crossover in super- 
conductivity, antiferromagnetism^i , or even the metallic 
transport'*''^, refiecting the substantial role played by the 
orbital-correlation in the multi-orbital electronic systems. 



II. THEORETICAL METHOD 

The vacancy ordered iron selenide AyFe2-xSe2 is a sys- 
tem with one-fifth of Fe being taken off and forming a 
characteristic \/5 x \/5 pattern (Fig. [H^a)). The pres- 
ence of the iron vacancy changes the space group from 
14/mmm (1)4/1) to 14/m (C4/1). As a result, the four-fold 
rotational symmetry is retained without the reflection 
symmetry with respect to the xz, yz and the diagonal 
planes. The absence of refiection symmetry implies that 
there is no symmetry axis. Hence there is an extra degree 
of freedom in defining the zero value of angles for basis 
functions of Cih- For instance, the basis function cos(46') 
becomes cos[4:{9 — 9o)] with 9o being an free parameter. 

To implement the I4/m symmetry for Aj^Fe2-xSe2, we 
have constructed a 12-orbital tight-binding model with 
four Fe atoms (labeled by A, B, C, and D) per cell and 
three t2g orbitals (d^z , dyz , and dxy with x and y referring 
to Fe-Fe directions) per Fe to investigate the magnetic 
instability'*^. The resulting FS is shown in Fig. [TJb) with 
two hole pockets, ai the one centering on (0,0) and Q!2 on 
(7r,7r), and four electron pockets, /3i the one and its in- 
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FIG. 1: (Color online) (a) The illustration of the Fe plane 
with the left chiral \/5 x \/5 vacancy ordering. The solid 
circles denote Fe atoms and the empty ones the vacancies. 
Every four nearest vacancies enclose the unit cells. Each unit 
cell contains four Fe atoms, labeled by A, B, C, and D. x and y 
are the primitive vectors, (b) The Fermi surface results from 
the 12-orbital tight-binding model: two hole pockets (ai and 
02) and four electron pockets (/3i and /?2 and their inverse 
partners). 216 fc-points on the Fermi surfaces are selected 
and shown by circles (zoom in to see). Points on FS-ai, 02, 
/?i and /32 are parameterized by the angles, dai, Oa2, 6 pi, and 
6f)2, respectively. 



version around ±(7r/2,7r/2) and p2 around ±(— 7r/2,7r/2). 
To calculate the SC gap, we shall divide the momentum 
k space into a lattice with 200 x 200 points and approx- 
imate states on the FS sheets by picking totally 216 k 
points with 44 points on FS-ai, 40 points on FS-q;2, and 
66 points both on FS-/3i and (32- These k points on FS 
sheets are characterized by angles, Oai, 9a2, 9 pi, and 9^2 
as illustrated in Fig. [T](b). All angles are relative to the 
horizontal kx axis except that 6^2 is relative to the verti- 
cal ky axis due to that FS-/32 is a 90° rotation from FS- 
Pi . On the average, the energy resolution in the k space 
with 200 X 200 points within a 50meV window relative to 
the Fermi energy is 0.6meV, while for the chosen 216 k 
points, the momentum resolution is about A/c = tt/IOO 
and we estimate the corresponding energy resolution is 
about 4meV. 

We shall assume that the interaction between electrons 
on each Fe atom is given by the generalized Hubbard 
model. 



Hi = \ UYfnai.itnai,ii 

i I L a=l 



a,b(a>b) 



U' ^ I naIAnbI,i — 2JH^aI,i ■ ^bl, 



+Jc {diI^,^d\^^^^dM,^idM,i^ + h.c^ | , 



(1) 



3 



where / is the index for four Fe atoms and a is the index 
for three orbitals. It was found''^ that as the interac- 
tion turns on, the spin fluctuation is strong at q —{tt,tt) 
and (tTjO), in which the former is related to the checker- 
board antiferromagnetism and the latter is related to the 
stripe antiferromagnetism. These two magnetic states 
competes each other and the preference between them is 
controlled by the Hund's rule coupling, Jh- 

The strong spin fluctuation is a possible pairing mech- 
anism. To investigate the superconductivity arising from 
the exchange of spin and charge fluctuations, we follow 
the FLEX approach^^ii^, in which the effective singlet 
pairing scattering matrix is given by 



r^^(k,k';c.) = 



k', wlFc 



(2) 



cd 



where w is a real frequency and a, h, c, and d as before are 
orbital indices. In Eq. the spin and charge vertices 
Ts and are F^'^^ = U, TY'^^''' = [/', F^-^'^' = 
Jh, TY'-"-'^ = Jc, and F^-^^ = U, F^''^^' = -U' + 
2Jh, VY'^'^' = 2[/' ~ Jh, F^'-^'^ = Jc, where non- 
vanishing vertices are only between the same Fe, and 
T labels orbitals and r 7^ t'. In the following, we will 
take the relations U' — U ~ 2Jh and Jc — Jh- The 
random phase approximation (RPA) for spin and charge 
susceptibilities are given by 

XRPA{f\,^) = xo(q,w)[i-F,xo(q,t^)]~\ (3) 
XRPAi^^^) = xo(q,'^) [1 + rcXo(q,w)]"^ , (4) 

with the bare susceptibility being given as 
1 



qML(k + q) 



(5) 



[x^{ci,^)C = }^ Ac,(k)A:^(k)^.(k- 

^ np [E^{\<i)]~nF [g.(k + q)] 
^ w + £;^(k)-£;^(k + q)-Ki(5' 

where is the orbital-band transformation matrix 

from V''a(k) = ^^a/i(k)7;i(k) and is the Fermi-Dirac 

distribution function. In this work, we shall set the tem- 
perature to zero. Here xo(q, is calculated by using 
the bare Green's function not the dressed one. The low 
energy physics of the scattering matrix is projected to 
scattering among the FS sheets. We shall denote the 
scattering matrix with a Cooper pair from (k',— k') on 
YS-v scattered to (k,-k) on FS-/i by r^^(k,k')ii, 



F^.(k, k') = <! ^ ^:^(k)A;^(-k) (6) 

xF^^(k,k';0)A.(k')Ac.(-k')}. 

Since we consider the even-parity pairing, we shall 
symmetrize the scattering matrix as F™'^"(k, k') = 

i [f^,(k,k')+f^.(k,-k')l. 
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FIG. 2; (Color online) The weighted scattering matrix 
[?;F(k)iiF(k')]~^ rji^''"(k, k') for k, k' being parameterized by 
Qai, da2, 6/31, and as defined in Fig. [T] Here the dotted 
white lines separate different FS sheets. Values of 9 range 
from 0° to 360° in each block bounded by dotted white lines, 
which represents the change of angles of k as it moves around 
the corresponding Fermi surfaces. Values of parameters are 
set with t/=1.0eV and Jh=0.2U. 



For a given gap function (7(k), a dimensionless coupling 
constant is defined by^ 
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E,,.Ek,k'^7Ik)WT5(k)rrr(k,k')5(k') 



k vf (k 



■<?2(k) 



(7) 

Here k and k' are restricted within FS fi and v, re- 
spectively. UF(k) = |Vk-B;j(k)| is the Fermi velocity, 
and Np (=216) is the number of k on Fermi surfaces 
we choose. The coupling constant includes contributions 
from pocket-pocket scattering processes, and it is helpful 
to extract them and to denote by a matrix (A)^!/, satis- 
fying A — J^fiuWt^i'^ where /x, i/ refer to FS sheets, /i, 
h'—l to ai, 2 to a2, 3 to and 4 to /S^. 

The stationary condition in Eq. (O (SX [g(k)] /(5g(k) = 
0) leads to the eigenvalue problem 



1 



F™;"(k,k').ga(k') = A,.9,(k), (8) 



where the subscript a stands for different solutions. For 
a given eigenvalue Aa, its eigenfunction (7a(k) determines 
the symmetry of the gap, which could be s-wave or 
(i-wave. An s-wave state is characterized by ga(k) = 
ga{1Zk.) and rf-wave by ga(k) = —gailZk), where TZ is 
the 90° -rotational operation on k. Among the s-wave 
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solutions, we pick up the largest one and set it to be As 
and its eigenfunction to be g^. Similarly, for the rf-wave, 
they are denoted by Ad and gd- The eigenvalue equation 
Eq. (O is not identical to the Bethe-Salpeter equation 
in which magnitude of unity for the eigenvalue stands for 
the formation of superconductivity. Instead, eigenvalues 
in Eq. (O stand for the pairing strength. Therefore, the 
SC state is determined by largest one of A^ and A^. 



III. NUMERICAL RESULTS 

A. Scattering Matrix 

The scattering matrix shows some clues to the gap 
function. In order to gain more condensation energy, 
the gap should be larger in the region where the scat- 
tering matrix is maximal in magnitude. As the interac- 
tion is repulsive, the gap becomes anisotropic and sign- 
changing. Fig. [2] shows the "weighted" scattering matrix 
[ui.(k)wF(k')]"^f™'="(k,k') for [/=1.0eVand Jh=0.2U. 
Here because the inverse Fermi velocity (IFV), l/vpCk), 
is related to the density of states (DOS), it is used as 
a weighting factor at each point. The scattering matrix 
contains intra and inter pocket scattering channels and 
different pockets are discriminated by the white dotted 
lines. The angles, 9^1, 0a2, 0pi and 9^2 parameterizing 
momenta k (k'), increase from 0° to 360° as the coordi- 
nate increases from the bottom (left) to the top (right). 
Although we only show the result from one of the electron 
pocket /3i due to the nature of being even parity for 
r™'^"(k, k'), the result at the inversion point can be de- 
duced. 

In Fig. [21 the intra-electron-pocket (/3i-/3i or 
scattering is much stronger than other channels and gives 
peaks at two momenta, « 48° and 225° as the an- 
tipodes (on the major axis) of the electron pockets. The 
strong intensity of the scattering matrix at these two 
momenta is expected because large IFV happens about 
the antipodes of the elliptic FS (l/wF(k) is shown in 
Fig. [6]). One of the reasons why intra-electron-pocket 
scattering is strong is due to the much larger DOS on 
the electron pockets. This also explains why the inter- 
electron-hole-pocket scattering is stronger than the intra- 
and the inter-hole-pocket scattering. Furthermore, the 
intra-electron-pocket scattering comes from exchange of 
fluctuations at vectors q «(0,0) and (7r,7r), and the latter 
is driven by the substantial checkerboard antiferromag- 
netic fluctuations. The strength of DOS, however, can 
not explain why the inter-electron-pocket (/3i-/32) scat- 
tering is much small because the stripe antiferromag- 
netic fluctuations at q =(7r,0) are not weak at all. A 
detailed analysis shows that the reduction of the inter- 
electron-pocket scattering is due to the difference in the 
orbital-band matrix element. In Fig. [3l we show the 
orbital weight |yla;i(k)|^ for every FS sheet . Here lines 
represent Fe-A/C and legends with circles and triangles 
represent Fe-B/D. Note that different curves exhibited 
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FIG. 3: (Color online) Orbital weight \Aa^,{'k)f around FS 
sheets, /i = ai (a), Q2 (b), /3i (c), and /32 (d). Here lines 
(solid, dashed, and dash-dotted) represent Fe-A/C and leg- 
ends (triangular, circular, and square) represent Fe-B/D. Dif- 
ferent curves are related by rotational symmetry (see text). 



in Fig. 3 are related by rotational symmetry^. For ex- 
ample, a 90° rotation about the center of the unit cell 
located at the position of Se changes Fe-A to Fe-B and 
dxz to dyz- As a result, as shown in Figs. [3Ka) or [3Kb), 
weight curve of orbital-(yz, B) is equal to that of orbital- 
{xz,A) shifted by a phase 7r/2. Similarly, since orbital- 
{xz, C) can be obtained by a 180° rotation of orbital- 
{xz,A), they have the same weight too. This is the 
relation within the Fermi surface of ai or a2- For re- 
lations between different Fermi surfaces, since the same 
rotation changes FS-/3i to FS-/32 and relevant orbits also 
change under the rotation, the orbital components on dif- 
ferent FS-/3 would be highly different. For instance, in 
Figs. [2fc) and[2Kd), the major component on FS-/3i is 
orbital- {xy, B) , while it is (xy. A) on FS-/32 and weights 
of other orbits are much smaller. Because vertices in 
the fluctuation exchange formula in Eq. ^ only couple 
orbitals on the same Fe, the inter-electron-pocket scat- 
tering will be heavily reduced. We have examined the 
case of neglecting the orbital-band matrix elements in 
Eq. © by setting ^^^(k) = l/\/T2 and find that the 
intcr-electron-pocket scattering becomes comparable to 
the intra-electron-pocket one. Therefore, the small inter- 
electron-pocket scattering is due to the small magnitude 
of the orbital-band matrix element. 



B. Phase Diagram of Superconductivity 

After a better understanding of the scattering matrix, 
we now show the results of eigenvalues from the eigen- 



5 



J..^O.OU J,,^0.2U J^^OAU 




U(eV) U(eV) U (eV) 

FIG. 4: (Color online) Eigenvalues or dimensionless coupling constants As (s-wave) and Ad (rf-wave) as functions of U for 
Jh=0, 0.2 1/, and 0.4(7 from the left panel to the right one. Solid lines are for As and dashed lines are for Aj. Insets are the 
difference of eigenvalues, AA = As — Ad, to exhibit the favorable state. The eigenvalues for Jh=0 are smaller, but they grow 
up fast and solutions become unstable at about U=1.0eV. 



value problem in Eq. ([5]). Fig. U displays the largest 
eigenvalues of s-wave (As) and d-wave (A^) as functions 
of U in three Hund's rule couplings, Jh=0, 0.2U and 
0.4 1/. Xg is represented by solid blue lines while Xd is 
represented by dashed red ones. As the interaction U is 
increased, quantum fluctuations enhances the SC pair- 
ing and then A's increase. Moreover, we observe that 
at different Jh, the dimensionless coupling constants for 
different symmetry grow in different speeds. As shown 
in the insets where the differences of two eigenvalues 
AA = Xs — Xd are plotted, the d-wave is always more 
stable against the s-wave (AA < 0) at Jh=0, while the 
preference is reverse (AA > 0) at Jh—OAU. Between 
them at Jij=0.2[/, the rf-wave is favored for [/<1.0eV 
but becomes unfavorable above that. 

The above results indicate that the Hund's rule cou- 
pling is a parameter of controlling the SC symmetry as 
well as the AFM order—. This is further supported by 
Fig. [SI in which AA versus Jh at different U is plot- 
ted. Except for the small phase space with U =1.0eV and 
Jh < 0.4C/, AA increases monotonously with Jh, indi- 
cating that the Hund's rule coupling favors the s-wave. 
However, AA is not a monotonic function of U. There 
exists a critical Jh, Jh,c ~ 0-2U, separating two states; 
below it, AA decreases as U increases, but above it, AA 
increases with U. Unless further symmetry is broken, s- 
wave and d-wave do not mix and the phase transition 
between them is first order. 
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FIG. 5: (Color online) The eigenvalue differences between the 
s -wave and the d-wave, AA, as functions Jh/ U in four cases 
of U: U=0.7eV (solid line), C7=0.8eV (dash line), C/=0.9eV 
(dotted line), and U=1.0eV (dash-dotted line). A critical Jh 
(~0.2[/) is observed as the transition point between s-wave 
and d-wave. 



tors. We shall first examine the case of s-wave solution 
in Fig. [6] and then the d-wave one in Fig. [T] 



C. Gap Functions 

With the phase diagram of superconductivity, one still 
needs to know the gap function explicitly to understand 
its physics, especially in this new type of system with 
space group different from other iron-based superconduc- 



Following Ref.[6,], we shall speak of s-wave by 
gaiTZk) — (?Q(k) with TZ being the 90°-rotational op- 
eration on k. Fig. [6] shows the s-wave gap function 
gs(k) along FS sheets ai, a-i, and for ?7=1.0eV and 
Jh=0-'2U. 5s(k) on FS-/32 not shown is the same as on 
FS-/3i when 6/32 = Opi. Blue squares are the eigenvalue- 
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FIG. 6: (Color online) s-wave eigenfunction (7s(k) for 
t/=1.0eV and Jh^0.2U. Panel (a) is along FS-ai, (b) along 
FS-a2, and (c) along FS-/3i. Squares represent the eigenvalue- 
problem solution and the solid lines are the fitting curves us- 
ing function forms shown in the text. Due to the s-wave 
character, the gap function on FS-/3i and FS-/32 is related by 
gs{0i32) = gs{Ofji) when 6^2 = dpi. The factor, l/uF(k), be- 
ing proportional to DOS, is also shown by the dashed green 
lines, a is the lattice length and its value from transmission 
electron microscopy^ is about 6.15A. 



problem solution, while the solid red lines are the fit 
guided by eye. One of the features is that such defined 
s-wave gap have multi- nodes on FS-ai, j3i and /32. More- 
over, because there is no reflection symmetry, there is an 
extra degree of freedom in defining origins of angles. As 
a result, the fitting curves for s-wave symmetry, chosen 
from lowest harmonic functions, are found to be 

g,{eai) = 0.031 COS [4 (6'c«i - 1°)] +0.0153, 

5s(^a2) = 0.028 cos [4(^,2+5°)] 

+ 0.008 cos [8 (6'„2 - 22°)] - 0.136, 

gs{epi) = 0.13[cos(fc;)-Kcos(fc;)] 



0.28 cos(fc;')cos(fc") -0.03, 



(9) 



y - sin(32°) 



where k'^ = k^ cos(32°) — ky sin(32 
A:yCOs(32°) and k'^ = k^coa{4°) + kySm{4°), k'^ 
— fca;Sin(4°) + fcyC0s(4°). Here cos40, cos80, cos(fca;) 



cos(fcy), and cos(fca;) cos(fcy) are the common s-wave 
bases. (Because FS-/3i is not located on symmetric point, 
a different type of bases is adopted.) However, due to 
broken reflection symmetry, origins of angles get shifted 
by 1°, -5°, 22° and 20°. These angles are the allowed de- 
grees of freedom in the C^h symmetry and their precise 
values are to minimize the total energy of the system. 

Similar to the s*-wave proposed in iron-pnictide 
superconductor—, where the SC order parameter on hole 
pockets and on electron pockets have opposite signs to 
make superconductivity stable, here we have the similar 
mechanism. In Fig. [51 we see that the mean value of 
the gap on FS-Q!2 has an opposite sign to those on FS-ai 
and FS-^'s, which suggests these inter-pocket scatter- 
ings gain energy. Due to their larger scattering strength, 
the a2 — /3's scattering processes gain most energy. The 
ai —/3's scattering processes are not favorable due to large 
cancellation by the oscillating behavior between positive 
and negative of gap on FS-ai. For the intra-pocket scat- 
tering processes, because the Fermi surfaces are small and 
if the fluctuations are smooth over the wave vector q, the 
interaction would not able to change from repulsive to at- 
tractive. However, the repulsive strength can be reduced 
by making the gap oscillatory and even with the sign 
being changed. The reason why gap functions take the 
particular forms of Eq. (jHl) can be traced back to the be- 
havior of the inverse Fermi velocity (IFV), which is shown 
as the dashed lines in Fig. [6] Clearly, the oscillatory be- 
havior of gap functions results from oscillatory IFV. On 
FS-/3i, the IFV is peaked at 6*^1 w225° where the gap is 
a peak too. The fact that the other point at Oj^i w48° 
with high IFV does not lead to a large gap is because the 
repulsion can be lowered in the intra- and inter-electron- 
pocket scattering processes. Due to that gap functions 
on FS-/? change slowly, it is clear that when the gap on 
FS-a2 is in-phase with IFV and dose not change sign, 
one gains energy in a2 — /3's scatterings but loses energy 
in the intra-pocket scatterings. The compromise between 
the intra- and inter-pocket scattering processes results in 
the gap function on FS-ai with the gap function being 
about ^ out of phase relative to IFV. In the following, we 

list numerical values of the coupling matrix Xg/Xg, which 
are consistent with the above explanations: 



-0.0275 0.2628 -0.1039 -0.1039 
-2.829 2 1.603 3 1.603 3 
-1.1398 -0.1947 
-1.1398 



where due to that is a symmetric matrix, the lower- 
triangular matrix elements are omitted for brevity, and 
As -E^.(A.W= 0.0837. 



2. d-wave 



Similar to the s-wave, the d-wave that we shall speak of 
obeys ga{TZk) = —ga(k) with TZ being the 90°-rotational 
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FIG. 7: (Color online) d-wave eigenfunction (;d(k) for 
i7=0.9eV and Jh=0.0. Due to the d-wave character, 
gd{0p2) = -gd{Opi) when 0^2 = fi'si- Here legends repre- 
sent the same meaning as those in Fig. [S] while the fitting 
curves are different and are discussed in the text. 



operation on k. In Fig. [71 we show the d-wave gap 
function gd(k.) for f/=0.9eV and Jh—0- 5d(k) on FS- 
/3i and on FS-/?2 is equal in magnitude but has opposite 
signs, i.e., gd{dp2) = ~gd{6f}i) when 6*^2 = dpi. It is seen 
that the gap function of d-wave always exhibits nodes on 
any Fermi surface. The fitted functions in Fig. [7] (solid 
red) are given by 

ffd(0ai) = 0.052 cos [2 12°)], 

Sd(0a2) = 0.18 COS [2 (0„2- 39°)], 

gdiOpi) = G.04[cos(fc;)-cos(A:;)] 

- G.3sin(2fc;')sin(2fc;'), (10) 

where k'^ = cos(20°) - ky sin(20°), k'y = k^ sin(20°) + 
A:ycos(20°), /c^' = A::^ cos(2°) + fcy sin(2°), and fc(; ^ 
— fc:rSin(2°) + fcj,cos(2°). Note that instead of cos46' 
(which is correct for the s-wave), d-wave is represented 
by cos 29. Similar to the s-wave case, although nodes ap- 
pear on every FS sheet, they are shifted away from the 
diagonal directions. 

The main cause for the occurrence of d-wave differs 
from that for the s-wave and is primarily due to reduc- 



coherent 
stripe AFM 
s-SC 



incoherent 
checkerboard AFIVl 
d-SC 



FIG. 8: Phase tendency of iron-vacancy-ordered AyFe2-a:Se2 
versus the Hund's rule coupling over the Hubbard U (Jh/U) 
and the pressure {P ). Here arrows point to the direction that 
the corresponding parameters increase their strength. Since 
increase of the pressure generally reduces Jh/U (see text), 
the increase of Jh/U or the decrease of P would drive the 
multi-orbital electronic system from coherent to incoherent, 
or from stripe to checkerboard AFM state, or from d-wave to 
s-wave SC state found in this work. 



tion of repulsion in intra-pocket and inter-electron-pocket 
scatterings through its sign-changing character . As a re- 
sult, both antipodes on FS-/3i {9fsi «48° and 225°) have 
large gaps in Fig. [71 in comparison to a single hump 
of the s-wave in Fig. [HI The coupling constant matrix 
(Ad)p^ summarizes the above effects 



Xd 
Xd 



0.0033 0.0361 0.0066 0.0066 
-0.0393 0.262 3 0.262 3 
-0.0459 -0.0066 
-0.0459 



where Ad = ^{Xd)fj,u — 0.0305 . Note that not only 
the repulsion between intra-pocket electrons is reduced, 
due to the presence of nodes, the attraction between 
inter-pocket electrons is also reduced. Although the d- 
wave state might not gain much energy from inter-pocket 
scattering processes relative to the s -wave state, it saves 
more energy from the reduction of intra-pocket repulsion. 



IV. SUMMARY 

Proximity of antiferromagnetism and superconductiv- 
ity in Aj^Fei.6Se2 implies that the SC state could be de- 
rived from spin fluctuations. Under this assumption and 
assuming that superconductivity exists in the \/5 x a/5 
iron- vacancy ordering state, we study the SC states from 
the effective pairing interaction in the FLEX approxi- 
mation based on our previous 12-orbital tight-binding 
model^. Similar to that of iron-pnictide superconduc- 
tors, for the spin-singlet superconductivity, s-wave and 
d-wave states are found to be close in energy. In par- 
ticular, a quantum critical point is found at Jh ~ 0.2U, 
below which the d-wave prevails, while above which the 
s -wave wins over. Furthermore, unlike the iron-pnictide 
superconductors, the iron- vacancy order lowers the sym- 
metry to the group Cih so that the gap functions do not 
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need to be reflection symmetric. Therefore, symmetric 
axes of the pairing momenta do not need to coincide with 
axes of the unit cell. As a result, the relative orientations 
of the pairing wave functions to the lattice, i.e., origins 
of angles for the pairing momenta, become new degrees 
of freedom for characterizing the superconducting gaps. 
This implies that the complete order parameters for char- 
acterizing the superconductivity in Aj,Fei.6Se2 consist of 
both the magnitude of the gap and the orientation of the 
pairing wave function relative to the underlying lattice. 

Finally, while our work focuses on the superconduc- 
tivity, the finding of quantum critical point controlled 
by Jh in superconductivity is not an accident. In the 
magnetism, we have shown that Jh controls the mag- 
netic phase transition from the stripe to the checker- 
board antiferromagnetism"*^. In addition, Haule and 
Kotliar— show that the crossover between coherence 
and incoherence is also determined by Jr- Therefore, 
the Hund's rule coupling plays the dominant role in 
determining localization, magnetism, and superconduc- 
tivity of iron-based superconductors. Experimentally, 
it is observed that applying pressure drives the mag- 
netic phase transition^SiSi that accompanies the metal- 



to-semiconductor transition. Although in general, apply- 
ing pressure affects itinerary electrons and could change 
Jh and U, the experimental results, when combined 
with our results on magnetism^, indicate that increase 
of the pressure generally reduces Jh/U in this system. 
Therefore, the phase tendency for iron-vacancy-ordered 
AyFe2_a;Se2 can be summarized in Fig|8]. The tendency 
implies that by appropriately changing the applied pres- 
sure, transition between different pairing symmetry of 
superconductivity can be induced. While the phase ten- 
dency shown in Fig. [8]needs further supports from exper- 
imental confirmation, our results provide a useful guide in 
clarifying the origin of unconventional superconducting 
gaps and searching for a possible quantum phase transi- 
tion in iron- vacancy-ordered Aj,Fe2_a;Se2. 
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